Particle motion during the compaction of granular matter 
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We track particle motions in a granular material subjected to compaction using a laser scattering 
based imaging method where compaction is achieved through thermal cycling. Particle displace- 
ments in this jammed fluid correlate strongly with rearrangments of the Voronoi cells defining the 
local spatial partitioning about the particles, similar to previous observations of Rahman on cooled 
liquids. Our observations provide further evidence of commonalities between particle dynamics in 
granular matter close to jamming and supercooled liquids. 



PACS numbers: 45.70.Cc,81.05.Rm,45.70.-n 

Granular materials are strongly interacting particle 
systems that admit an admixture of collective particle 
motions, as in crystals, and random particle motions, 
as in simple fluids. Under flow or deformation condi- 
tions these materials pass intermittently between a solid 
"jammed" state and a state of flow. When a granular 
material jams, the individual particles are in stable me- 
chanical equilibrium with their local neighbors [1 . Small 
perturbations such as tapping [5j or shearing [3j can move 
the particles, and lead to the evolution from one jammed 
state to another, generally more compact state. The in- 
crease in particle density with the number of perturba- 
tions is similar to the slow thermally driven dynamics 
of fluids close to their glass transition. Observations of 
phenomena of this kind has led to the development of a 
general framework for the jamming transition, with tem- 
perature, density, and external forcing as the three main 
variables characterizing this dynamical transition [TJS]. 

Although there have been a number of recent studies of 
particle trajectories in glass- forming colloidal fluids 
and other strongly interacting particle systems, little is 
known about the nature of particle trajectories in gran- 
ular materials. While recent studies have successfully 
observed particle positions in a 3D jammed granular ma- 
terial [71 |S], particle trajectories were not observed. It 
seems plausible that the particle motion has features in 
common with glass- forming liquids, such as the emer- 
gence of string- like collective motion [9], and based on 
this hypothesis we examine whether the shape fluctua- 
tions in the Voronoi cells environments about the par- 
ticles in compactified granular media correlate with the 
direction of particle displacements, as noticed by Rah- 
man [TO] in his pioneering studies of collective particle 
motions in cooled liquids. On the other hand, the non- 
equilibrium nature of these driven particle systems leads 



us to expect differences from thermalized liquids that re- 
quire a better understanding. In this paper, we provide 
the first direct measurements of 3D particle trajectories 
during compaction, with a focus on the statistics of par- 
ticle motion. 

In particular, we study particle rearrangements in a 
column of granular materials during the course of ther- 
mal cycling. We use a laser sheet scanning method to 
reconstruct the configurations of the system at the end of 
each cycle. The motion of individual particles is tracked 
over 10 cycles, and we analyze how the local geometry 
affects local particle dynamics. 

Thermal cycling provides a method of granular com- 
paction without exciting strong motions of particles. By 
thermal cycling we meanalternately heating and cooling 
the material and container, with one heating and cooling 
cycle corresponding to one perturbation. This technique 
leads to exponential compaction of the material with two 
characteristic timescales if the granular material and 
container have different thermal expansion coefficients. 
While most methods of forcing rearrangements of gran- 
ular matter, such as tapping of the boundaries or air 
fluidization agitate individual particles significantly, i.e. 
provide particles kinetic energy relative to their neigh- 
bors, thermal cycling provides a quasi-stationary method 
of evolving the structure towards a jammed state. 

Experimental Setup: The granular material used is 
soda-lime glass beads (diameter D = 3 mm ± 0.3 mm), 
poured into a transparent polymcthylpcntcnc (PMP) 
cylinder of diameter 5cm (16.6 D) to a height of 19.5 cm± 
0.5 cm (63.3 D to 66.6 D). This is the same materials as 
used in previous work on thermal cycling [11 . The ther- 
mal expansion coefficient for PMP is 1.17 x 10" 4 K" 1 , 
about one order of magnitude larger than that of the 
beads (9 x 10~ 6 K _1 ) [TT]. The beads were immersed in 
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index-matching oil (Cargille Labs Type DF [20]) which 
contains 2.08 pg/mL laser dye (Nile Blue 690 perchlo- 
rate). A transparent weight of 145.5 g was placed on top 
of the particles to monitor the filling height and to apply 
a controlled vertical force. 

The cylinder is placed inside a transparent water tank 
for thermal cycling at a rate of 16° C per hour. As 
the PMP container was heated/cooled the cross-sectional 
area expanded/contracted by 0.46 mm 2 for every degree 
heated/cooled (6.5 % of a single particle cross-section). 

While bead expansion was negligible, the sample and 
oil are index-matched at room temperature. Since the 
index of refraction changes with temperature differently 
for glass and oil, we can only image our sample at the 
end of each thermal cycle. 




FIG. 1: Experimental setup: Glass beads immersed in index 
matching oil in a cylinder are subjected to thermal cycling 
via a water bath. The fluorescently dyed fluid is imaged via 
a laser sheet that is moved across the sample. One of 600 
cross-sectional images shown. 

Three dimensional images of the system at each cycle 
were obtained using a laser sheet scanning method, first 
described in [12]. Briefly, a laser sheet is sent through 
the specimen at the excitation frequency of the laser dye 
to produce a cross-sectional image of the granular mate- 
rial in which the particles appear as dark circles within 
the dyed fluid (see Fig. [IJ. Images are taken with a high 
sensitivity cooled CCD camera (Sensicam, PCO). To ob- 
tain a three dimensional image, the laser is translated in 
100/im increments along the axis of the cylinder for a dis- 
tance of 60 mm (20 D), starting from the top of the pile. 
The camera is not translated, but the changing object 
distance is accounted for in the analysis. 

A 3D bandpass filter is used to reduce noise and to 
smooth out images of the beads. A particle tracking 
algorithm by Crocker [T3] and Weeks [TJ] detects the 
centers of > 98 % of all beads to within 100 /xm, which 



is smaller than the uncertainty in the bead diameters. 
The remaining < 2 % of bead positions are detected 
manually with comparable resolution using an interactive 
3D visualization of images with an overlay of extracted 
particles. This allows us to generate complete maps of 
3D jammed states. Particle rearrangmenets are followed 
through multiple thermal cycles. The particles move by 
less than 0.5 D in each thermal cycle, a scale small enough 
to allow for reliable particle tracking. 
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FIG. 2: Packing fraction vs. cycle number for various tem- 
perature differentials 

Results: We apply ten thermal cycles for three differ- 
ent temperature differentials: 25° C, 40° C, and 60° C. 
The average packing fraction is determined from the ver- 
tical position of the flat glass weight on top of the pile 
in the 3D image. The initial conditions are generated 
by letting particles settle into the oil filled cylinder. Af- 
ter the particles settle, the glass top is allowed to sink 
onto the surface of the fluid immersed particles. Fill- 
ing in oil first is necessary to eliminate air bubbles from 
the sample, but leads to significant variability in the 
initial packing fraction. Nevertheless, the increase in 
packing fraction during thermal cycling can be fit with 
an exponential p(t) = po — A e'/ r as shown in Fig. [5] 
with po = 0.619 ± 0.0008, A = 0.028 ± 0.0007, and 
r = 4.87 cycles ± 0.31 cycles. The compaction rate is 
slower than the faster of the two time scales, 2.72 cycles, 
reported by |llj , consistent with a weaker forcing of the 
particles due to buoyancy forces. 

We find that the pair correlation function g(r) does 
not evolve significantly during compaction, as shown in 
Fig. [3] This is not surprising since the particles only need 
to move closer by about 0.001 D on average to achieve an 
increase in packing fraction of 1.5%. The overall shape of 
g(r) is consistent with a random arrangement of particles 
in the presence of some ordering at the boundaries. 

The Voronoi construction defines a unique partitioning 
of space where each particle is assigned a local neighbor- 
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hood, the Voronoi cell, corresponding to the region of 
space closest to a given particle. The Voronoi cell struc- 
tures were determined using the QHULL algorithm |15j . 
To avoid boundary artifacts, Voronoi cells were computed 
for all particles, but only the volumes of particles that are 
not on the boundaries we used in the analysis. The dis- 
tributions of Voronoi volumes is shown in Fig. [4] The 
distribution is non-gaussian with a sharp cutoff at small 
volumes and a broad range of large volumes, similar to 
distributions observed in simulations in glass- forming liq- 
uids |16j and in recent experiments of granular materials 
8 . One possible functional form of the distribution of 
voronoi volumes proposed by Aste et al. [8] contains only 
one fitting parameter k. Fitting our Voronoi cell volume 
distributions to this function yields k = 24, apprecia- 
bly larger than the k value, k = 12, suggested by Aste 
et al. [8 , however. The larger polydispersity in bead 
diameters in the experiments presented here (10 % as 
opposed to 2 % to 5 %), could explain the wider distribu- 
tion of Voronoi volumes. We find that the distribution 
of Voronoi volumes does not change significantly as the 
system is compacted. 

Following Rahman's investigation of cooled liquids, we 
analyze particle motion and its relation to the Voronoi 
cell shape. We define the vector Vij to be the displace- 
ment of the center of particle i from cycle j to j + 1 (see 
Fig. [5]). Next, we define a vector describing the Voronoi 
cell shape relative to the particle center. From the per- 
spective of the particle, the vertices of the vornonoi con- 
struction correspond to the direction in which a particle 
"sees" a void between three neighbors. The mean po- 
sition of all vertices therefore indicates where a particle 
would sec more local void spaces. The vector w^j defines 
the displacement from the center of particle i at cycle j 
to the mean position of its corresponding Voronoi ver- 
tices, as shown in Fig. [5] The averages for \u\ and \v\ are 
approximately 1 mm (0.33 D) larger than the accuracy 




FIG. 5: Definition of particle displacement vector v and u 
vector for a cross-section of the 3D Voronoi volume. Since 
this is a cross-section of the 3D Voronoi cell, the particle need 
not be in contact with any of the edges of the Voronoi volume. 

with which we detect particle positions. 

To determine the correlation between particle displace- 
ment and Voronoi shapes, we then consider the alignment 
of Uij and Vij. Specifically, we calculate cos(#) for each 
image, where 9 is the angle between j and Vij for each 
particle. The probability distribution of cos 6 averaged 
over all cycles is shown in Fig. [6j In addition, we show 
the distribution for the first and the last cycle, and a 
control for two 3-D images without thermal cycling be- 
tween the images. We find that the distribution peaks at 
cos# = 1, indicating that the particle moves toward the 
center of the vertices of its Voronoi cell. Particle motion 
and Voronoi shape are indeed correlated. 

Discussion: 

By using the laser sheet scanning method, we visu- 
alize the internal structure and dynamics of a jammed 
granular material. Under quasi-static forcing via ther- 
mal cycling flows are slow enough (less than w 0.5 D 
between frames), so we can observe the microscopic dy- 
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namics of this strongly interacting particle system. While 
the differences in local structures are too small to be ob- 
servable in local distributions such as g(r) or the Voronoi 
volume distribution, the Voronoi volume appears to hold 
important clues for the future dynamical evolution of the 
particles. Voronoi reconstruction indicates that particles 
tend to move toward the centers of Voronoi vertices as 
long ago observed in Rahman's pioneering simulations of 
particle motion in cooled liquids [10]. This correlated 
displacement in our granular fluid is independent of the 
size of the Voronoi volume, the velocity of particles, and 
the position of particles within the granular column. 

Rahman [TO] and others following him [T7] inferred 
from the correlations in the tracer particle displacement 
and the Voronoi cell shape that the particles were moving 
in tube-like environments, thus anticipating the string- 
like motion in cooled liquids that has recently become 
widely appreciated [OJ [TB]- This raises questions about 
whether such collective interparticlc displacements also 
characterize particle motion in compactified granular ma- 
terials. The absence of equilibrium in granular fluids, 
however, complicates the definition of "mobile" and col- 
lectively moving particles in granular materials, since 
there is no model such as Brownian motion to provide 
a "baseline" for this comparison. A recent study of a 
quasi-2D air-driven granular material |19j indicates the 
presence of collective string-like motion as in cooled liq- 
uids so the general presence of such collective motions 
in driven granular media seem highly plausible. How- 
ever, we expect important differences to emerge in sys- 
tems such as ours where the volume of the material is 




FIG. 6: a) Preferred Direction PDFs for system with AT = 
40° C for each thermal cycle, b) Preferred Direction PDFs 
for system with AT = 25° C for each thermal cycle, c) Pre- 
ferred Direction PDFs for system at constant temperature d) 
3D reconstruction of system from the perspective of a given 
particle facing five vertices of its Voronoi cell. Notice that the 
vertex is located in the void between particles. 



evolving in time with compaction. We plan to study the 
evolutionary or aging dynamics of our granular fluid in 
the future. 
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